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Abstract 

We present the first calculations of proton elastic and inelastic scattering in 
which the Pauli blocking operator contains the leading non-spherical components 
as well as the usual spherical (angle-averaged) part. We develop a formalism for 
including the contributions to the effective nucleon-nucleon interaction from the 
resulting new G-matrix elements that extend the usual two-nucleon spin structure 
and may not conserve angular momentum. We explore the consequences of parity 
conservation, time reversal invariance, and nucleon-nucleon antisymmetrization 
for the new effective interaction. Changes to the calculated cross section and spin 
observables are small in the energy range from 100 to 200 MeV. 
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1 Introduction 



Proton elastic and inelastic scattering at energies above about 100 MeV are usually 
described by distorted-wave calculations based on an effective nucleon-nucleon (NN) 
interaction. When the interaction with the projectile proton is summed, or "folded," 
over all the nucleons in the target, the resulting potential can be used to model elastic 
scattering. In addition, the same effective NN interaction becomes the transition po- 
tential in the Distorted Wave Impulse Approximation (DWIA) that connects to excited 
states of the target when the struck nucleon moves to a new shell-model orbit, creating 
a particle-hole pair. The many-body effects of the nuclear medium are usually incor- 
porated through modifications to this effective NN interaction that depend on the local 
nuclear density. 

Systematic studies [TJ EH HI E| have shown that one important contribution to the 
many-body effects is Pauli blocking, particularly at the lower end of the intermediate 
energy range. This is the mechanism which prevents nucleons in the nuclear medium 
from scattering to occupied intermediate states jB] • This restriction is included through 
a projection operator in the Bethe-Goldstone equation for the G-matrix elements that 
describe the effective NN interaction inside the nuclear medium. 

The usual practice is to average the Pauli projection operator over the intermediate 
state scattering angle. If, instead of this "spherical" approximation, the non-spherical 
components are retained, new G- matrix elements appear |B1 El EI! ■ While remaining 
diagonal in total spin S and isospin T, the angular dependence allows coupled G-matrix 
elements that connect partial wave states where J ^ J' and I ^ I' beyond the \£ — £'\ = 2 
coupling generated by the tensor interaction. The expanded set of G-matrix elements 
also depends on M, the magnetic quantum number of the total angular momentum J. In 
studies that considered the effects on nuclear binding energies, small but non-negligible 
changes were found when these modifications to the G- matrix were included [HI El- It has 
also been suggested that the spherical approximation is adequate for the central and 
spin-orbit parts of the effective interaction [7j, those pieces that matter the most for 
elastic scattering and the excitation of natural-parity transitions. But no calculations 
have been made to check quantitatively in a comparison with scattering measurements 
how important it is to treat the non-spherical Pauli blocking components. 

In a recent calculation of the non-spherical G-matrix ^U] , we observed that the new 
couplings in J and i generated only small matrix elements, but the variation of the 
G-matrix elements with respect to M, the projection of the total angular momentum 
J, was comparable to the typical size of conventional medium effects. Thus it seemed 
appropriate to investigate whether these modifications could have an impact on nuclear 
reactions, and in particular on their spin observables, which are most sensitive to the 
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non-spherical components of the nuclear force. Proton scattering offers a rich set of 
polarization observables, including polarization transfer, that reflect the spin dependence 
of the effective nucleon-nucleon interaction itself. 

The dependence on M has been considered previously in calculations of deuteron 
binding when the deuteron is treated as a projectile traveling through the nuclear 
medium ^T]. While Pauli blocking reduces the deuteron binding energy as the den- 
sity increases, the amount depends on whether the projection of the deuteron's spin, 
\M\, is or 1 when the quantization axis is taken to be along the direction of motion of 
the deuteron through nuclear matter. This difference generates a Tp-type momentum- 
dependent tensor potential in the deuteron-nucleus optical model. Searches for such a 
potential have not been definitive because of interpretive complications arising from the 
strong coupling to breakup states [T2] . 

In this paper, we begin with the G-matrix described in Ref. ^U]. The important 
features of the calculation of the G-matrix elements are reviewed in Section II. In order 
to use this effective interaction in distorted-wave calculations with presently available 
computer programs, we must transform the G-matrix elements to a coordinate-space rep- 
resentation using a sum of Yukawa functions. In Ref. (3J these coefficients in the Yukawa 
expansion are fit directly to the values of the G-matrix elements. This is equivalent to 
earlier methods in which the nucleon-nucleon scattering amplitudes were calculated as 
a function of scattering angle or momentum transfer and then the Yukawa expansion 
ranges and coefficients were chosen to best reproduce these angular distributions [T3]. 
For technical reasons, we will follow the second scheme here (see Section III). In the 
process, we will introduce a multipole expansion of the new G-matrix elements in which 
the lowest order recovers the result for a spherical Pauli blocking operator. For the 
spherical Pauli operator, a comparison of the two methods shows practically identical 
answers. 

In Section IV we will compare the non-spherical and spherical treatments of the 
Pauli blocking operator for representative nuclear transitions at 100 and 200 MeV. We 
will include a comparison to the free, or density-independent, effective interaction and 
to an effective interaction that also contains relativistic effects. These comparisons will 
help gauge the importance of the non-spherical treatment of Pauli blocking relative to 
density- dependent effects more typically included. We will also include measurements of 
the cross section and analyzing power to illustrate these effects in relation to the quality 
of the reproduction of the data. Pauli blocking has its largest effects on the isoscalar 
central and spin-orbit terms in the effective interaction. These terms are well tested by 
a comparison to elastic proton scattering or transitions to natural-parity excited states. 

Pauli blocking is only one process that is important in the calculation of the effective 
interaction in the nuclear medium. Others, such as the effects of strong relativistic mean 
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fields [TUJ O CH] and coupling to A- resonances (TB] , increase the repulsion in the nuclear 
medium just as does Pauli blocking [T7j. In a complete treatment, these should be 
properly considered, along with the attraction expected to arise from many-body forces. 
Thus a set of calculations based on the conventional Brueckner-Hartree-Fock approach 
to nuclear matter only (such as those we present here), should not be expected to provide 
the final answer. However, a critical evaluation of any of these medium effects requires 
that the treatment of Pauli blocking not introduce systematic errors large enough to 
affect our interpretation when agreement with data is considered. We will show that at 
intermediate energies the inclusion of non-spherical components in the blocking operator 
produces changes that are modest in size compared to the main density-dependent effects 
typically included in the effective interaction. 



2 Calculation of the G-matrix elements 

The Brueckner-Bethe-Goldstone equation [THl EH I2DI EE] describes the scattering of two 
nucleons in nuclear matter. The presence of the (infinite) nuclear medium is included 
through Pauli blocking and a mean field arising from the interactions with all the other 
nucleons. It is convenient to express the momenta of the two nucleons, ki for the 
projectile and k2 for the struck nucleon, in terms of the relative and center-of-mass 
motion as 



k = (k x -k 2 )/2 

P = (k x + k 2 ) . U 

The total or center-of-mass momentum P is conserved in the scattering process. 

In analogy with free-space scattering, the nuclear matter Bethe-Goldstone equation 
is given by 

G(k<, k, P, E 0) = VP, k) + / ^ViK k'O ^^ ^ k--, k, P, E ) , (2) 

where V is the two-body potential. The energy of the two-particle system, E (with 
Eq its initial value), includes kinetic energy and the potential energy generated by the 
mean field. The latter is determined in a separate self-consistent calculation of nuclear 
matter properties and is conveniently parametrized in terms of effective masses [5]. At 
this point the direction and magnitude of P have not been specified. 

The Pauli projection operator Q selects intermediate states where both nucleon 
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momenta lie above the Fermi momentum kp: 



Q(k,P,A*) = {i k F (3) 

1 otherwise 

Visualizing the (sharp) Fermi surface as a sphere of radius hp, the condition above 
imposes the requirement that the tips of the k x and k 2 vectors lie outside the sphere. 
For applications to real nuclei, hp is treated as a function of the local nuclear density 

P(r) = ^ • (4) 

The matrix elements for the non-spherical Pauli operator may be written in a partial 
wave basis as 



( (£'S) J'M'\Q(k,P,k F )\{£S)JM) 

(£'rn e ,Sms\J'M')(JM\£m e Sms}(£'m e ,\Q(k,P,kF)\£m e ) , ( 5 ) 

where 

(e'm e/ \Q(KP,k F )\£m e ) = J dtl ^(fi) Y^(O) 0(|ki| - k F )<d(\k 2 \ - k F ) . (6) 

The step functions destroy the orthogonality that would otherwise exist for the spherical 
harmonics in the integral. This allows couplings where £ ^ £' and, through the recoupling 
coefficients in the summation of Eq. (5), couplings where J ^ J'. 

If the quantization axis for the projection quantum numbers in Eq. (6) is chosen to 
lie along P, then there is no dependence on azimuthal direction in the functions, and 
integration over dfl gives me = m^>. Because Q(\t,P , kp) is diagonal in ms, Eq. (5) 
gives M = M' . It then follows from Eq. (2) that G is diagonal in M. 

Eq. (5) can be further reduced using standard angular momentum algebra to give 

( (£'S) J'M | Q (k, P , k F ) | {£S) JM) 
= \Y.(- l Y +S+J J'Ll£'{£'^L l {]\m)W{J , J£ , £-L 1 S) (7) 

2 Li 

x (J'M, LiO\JM) 9 Ll (P,k,k F ) 

where 6L x {P,k,kp) is the combination of step functions defined in Eqs. (26)-(27) of 
Ref. [IT] . Here and elsewhere we use the notation X = \J2X + 1. In Eq. (7) the 
quantum number M is the projection of J along the direction of P. We have checked 
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that values of selected matrix elements calculated with this formula and Eq. (6) agree 
satisfactorily. 

The G-matrix elements are obtained by solving the integral Bethe-Goldstone equa- 
tion, Eq. (2), in the £, S, J, M, T basis. For any choice of P, the G-matrix elements, 
(£'J'M'\G ST (P)\£JM), may be expanded in spherical harmonics Y^a(P) where L is the 
angular momentum that recouples J to J'. The coefficients in the recoupling expansion 
are G LTS (£' J' JJ, P). Thus, 

(£'J'M'\G ST {P)\£JM) = V^^2(J'M',LA\JM) L Y LA (P) G LTS {£'J',£J,P) , (8) 

LA 

where the arguments of spherical harmonics are always unit vectors in the direction of 
the vector indicated. In the limit of a spherically averaged Pauli blocking operator only 
the L = terms in Eq. (8) survive. Equation (8) can be inverted to yield the expansion 
coefficients 

G LTS (£'J', £J, P) = i ^2(J'M, L0\JM) (£'J'M:P\G ST \£JM:P) . (9) 

J 2 M 

where the relation M = M' has been applied. The G-matrix elements on the right hand 
side are defined with respect to a basis \£JM:P) and M refers to the direction of P. 
For the nuclear matter calculations described here the magnitude of P has been chosen 
to be k\. 

Note that if the matrix elements (£' J'M : P\G ST \£JM : P) are independent of M 
and diagonal in J (spherically averaged Pauli blocking) the right hand side of Eq. (9) 
automatically vanishes unless L = because of a property of the C-G coefficients: 

J2(JM,L0\JM) = J 2 5 L0 . 

M 

The anti-symmetrized G-matrix elements are obtained by subtracting from Eq. (9) 
the same terms but with an additional phase of (— y+ s+T . This projects out the set 
of matrix elements with £ + S + T even, leaving the set that normally describes NN 
scattering. Due to parity conservation, L takes on only even values. The first step 
in our transformation to coordinate space was to obtain these L-dependent expansion 
coefficients. This expansion is expected to converge quickly in L, and the L = part 
was checked for consistency with the result for the spherically averaged Pauli operator. 

The G-matrix elements that we will use here were generated from a free-space NN 
interaction that is a modified version of the Bonn-B potential ^7] . Details can be found 
in Ref. [Sj- The model parameters were adjusted so as to achieve a good reproduction 
of the phase shift analysis results from the Nijmegen group [22] up to 325 MeV. 
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The number of matrix elements (£'J'M:P\G ST \£JM:P) that are coupled in Eq. (2) 
increases with M as it becomes possible to incorporate larger values of J. However, the 
M-dependence decreases with increasing J and we were able to ignore these effects 
on partial waves with J, J' > 6, which were calculated using the usual angle- average 
approximation. 

In general, the new elements of the G-matrix will introduce spin operators that 
depend on the direction of P, the sum of the projectile momentum k x and the momentum 
k2 of the particle encountered in the nuclear medium. For practical reasons it is desirable 
that P be fixed for a particular transition and incident energy. At the intermediate 
energies of interest here, when the incident laboratory momentum can be considered 
large compared with the momenta of the target nucleons, a natural approximation for 
P is kx- We shall show however that this is not a good choice here because it leads to a 
G-matrix which cannot be expressed in terms of the standard set of spin operators (for 
example, <j\ • 02, {v\ + a 2) ■ h, Su((j), and SuiQ), where the momentum transfer q and 
the normal to the scattering plane n establish a coordinate system with Q = q x n). In 
this paper we choose P to be parallel to Q and have magnitude fci. This does lead to a 
G-matrix with the standard spin structure. Our choice for the magnitude is a standard 
one in calculations of medium effects. 

3 Transformation of the G-matrix to coordinate 
space 

A calculation of proton elastic or inelastic scattering observables using presently available 
computer programs requires that the information on the NN interaction in the medium 
carried by the G-matrix elements be converted into the amplitudes of the effective NN 
interaction. Usually, the amplitudes for each spin and isospin operator are expanded in 
a Yukawa series as a function of the momentum transfer. We will begin this section with 
a description of the general formulas that connect the expanded G-matrix of Eqs. (8) 
and (9) to NN scattering amplitudes, regardless of the complexity of the coupling. This 
makes them available for any possible future application. For the calculations we show 
here, we will restrict ourselves to the set of operators associated with free NN scattering 
since these are the conventionally used set. After the general formulas, we will show the 
forms actually used for our calculations. 

We first convert from the £, S, J, M, T basis to a representation characterised by 
definite incoming and outgoing momenta and intrinsic spin projections. The amplitudes 
in the two bases are related by 
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(k',Sa'\G T {P)\k,Sa) = (2tt) 3 £ i e ~ £ ' (£'m e >, So'\ J'M')(£m e , Sa\ JM) 



c 



x (J'M',LA\JM)Y e , mei (k')Y; me (k) (10) 
x LV^ Y la (P) G LTS (£' J' , £J, P) , 



with a summation that runs over ( = £' J' M' £J Mmi'iniLh. The momenta k and k' are 
nucleon momenta in the NN center-of-mass system [see Eq. (1)]. 

With the normalization factors as in Eq. (10), the partial wave matrix elements 
G LST (£'J, £J, P) satisfy unitarity relations which, in the case of uncoupled partial waves 
in free space, have the form 

ih)-& ™ 

where \x is the NN reduced mass. This relation is equivalent to 

G e = - 8?r r exp(i5 e sin 5 e ) (12) 

where 5e is a real NN phase shift. This is the same definition of the partial wave G- matrix 
elements that was used in Appendix A of Ref . jH] . 

We denote antisymmetrized matrix elements by G T where 

(k', Sa'|G T (P)|k, So) = 2 (k', 5a'|G T (P)|k, So) . (13) 

For this to be correct, the summation of Eq. (10) and subsequent summations in this 
paper must restrict I and £' to satisfy 

For the purpose of applications in DWBA calculations we wish to separate the parts 
of G associated with the central, spin-orbit, and tensor operators usually used to describe 
the spin structure of the on-shell (k = k') NN scattering amplitude. In particular, we 
want to consider the form 

(k', So'\G T (P)\k, So) = ]T(k', So'\G LT (P)\k, So), 

L 

where 

(k' J( S<r'|G iT (P)|k, So) = G^ TS (6) 

+S S i[G L L f(e)(o l + o 2 )-n + G^(6)S 12 (q) + G^(6)S 12 (Q)} , (15) 
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— * — * — * — * 

and where q = k' x — hi = k' — k is the momentum transfer, n is the normal to the 
scattering plane, and Q = q x n. Each of the Gf TS is a function of the scattering angle 
0. The subscript indicates the spin operator in the NN amplitude corresponding to that 
coefficient, using C for central (both S = and 5=1), LS for spin-orbit, and TD and 
TX for the "direct" and "exchange" parts of the tensor interaction. 

Quite generally, the G-matrix (k', Sa'\G LT (P)\k, Sa) can be expanded in terms of 
the complete set of spin tensors T ksqs (S) whose matrix elements are 

(Sa'\r ksqs \Sa) = k s (Sa, k s qs\Sa'), (16) 

where ks runs from to 2S. The coefficients in this expansion are G k J qs (k, k') and are 
given by 

Gg£(k, k', P) = Trace «k', ^'|G LT (P)|k, Sa) r ksqs ) . (17) 
These coefficients G^(k,k',P) are defined and normalized so that 

(k',Sa'\G LT (P)\k,Sa) = S- 2 £ (Sa'\rt sqs \Sa) G^(k,k',P). (18) 

ksqs 

We emphasize that for a general P the amplitude (k', 5V|£r (P)|k, Sa) of Eq. (13) 
does not have the form of Eq. (15) but will contain other terms, e.g., S^P). However, it 
is shown in Appendix A that if P is chosen so that it becomes P' = ±P under a rotation 
of 7r about (k+k'), then the form Eq. (15) is sufficiently general. The following formulas 
will assume that this is the case. The coefficients in Eq. (15) are then given in terms of 
the GfJ s by 

G L J S (6) = S- 2 G^(k,k',P) (19) 

G L £{9) = f E^(«)0- k ' P ) (2°) 

qs 

G L T 1 J(e)- 1 -G^(9) = ^£r4(q)G^,k',P) (21) 

~G^(e) + G L T m = ^EF4(Q)GgJ(k,k',P). (22) 

These results follow from the orthogonality of the T ksqs with respect to the trace oper- 
ation and the formulas (S = 1) 

Trace(r 1(?s S-A) = VSn Y lqs (A) (23) 

8^ 



/87T 

Trace(r 2g5 (3(S-A) 2 -2)) = 3d — Y 2qs (A) , (24) 
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where A is an arbitrary unit vector. 

Using the partial wave expansion given in Eq. (10) the trace in Eq. (17) can be 
calculated explicitly to give, 



r £ e h) 

G L k ™(k,k',P) = ^ e ~ e 'SJ 2 J'kL(-l) J - J '~ e+k ^J J' L 

C [S S k s 

x (keqi,LA\k s qs) (£'m if , £m e \k e q e ) 

x Y e , mef (k>) Y eme (k) Y LA (P) G LTS (£'J',£J,P) 



(25) 



where the sum runs over ( = £'J'£Jm£'mekeqeA with the restrictions for antisymmetriza- 
tion noted above. 

For our present application, we choose to express the components of these tensors 
in a right-handed coordinate system with z along k and y along n. We call this the 
'standard' coordinate system. We choose P to be parallel to (k+k') so that the restricted 
form of Eq. (15) is valid. In the limit that the reaction Q- value is small compared to 
the bombarding energy, the direction of P is toward 6/2. We also neglect terms with 
J' 7^ J because these matrix elements turn out to be small. Eq. (25) reduces to 



r £ e h 

Gg2(k,k',p||k+k') = j2^ e, sj%L(-^r e+kt \j J L 

c (s S k s 

x {k e qe, LA\k s q S ) (£' qi, £0\hqe) 

x Y Vqi {9, 0)£ Y LA (6/2, 0) G LTS (£'J, U, P) 



(26) 



The values of G LTS (£' J, £J, P) were calculated using Eq. (9) in the form 

G LTS (£J,£J,P) = i {(J0,L0\J0) (i'JM = 0\G ST \£JM = 0) 
J 2 1 

+2 ]T (JM,L0\JM) (£'JM\G ST \£JM)} . ^ 

M>0 

Eq. (26) was split into two summations for the sake of faster computation as 
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G^(k,k',P||k + k') = L^^'JH-^S G LTS (e'J,£J,P) 



X 



r £ £' ki ' 

x (k e q e , LA\k s q s ) (l'q e ,lO\k e qt) 
x 1^(0,0)^(0/2,0). 



(28) 



When expressed in terms of the components of the tensors G^J^(k, k', P| |k + k') in the 
'standard' coordinate system, Eqs. (19)-(22) reduce to 



PyLTS 



(0) 



G L L f{9) 



q-2 P^LTS 
° ^OO 



P^LVT 



TD (#) 
X (#) 



PiLIT 



6 
V2 
72 
V2 
72 



3(1 - cos 9)G^ T + 2v / 6sin^ G£ T + ^(3 + cos0)<5 
3(1 + cos 9)G^ T - 2^6 sin G^ 1T + ^(3 - cos 6)G 



LIT 
22 



LIT 
22 



(29) 
(30) 

(31) 

, (32) 



where the arguments of the tensor components on the right-hand side have been 
dropped for simplicity. 

An additional transform was needed at the end to replace states where S = 0, 1 and 
T = 0, 1 with the singlet and triplet spin and isospin operators customarily used by the 
distorted-wave programs. 

In Ref. |TUj, it was noted that, if the M-dependent G-matrix elements were simply 
averaged over M, the result was very close to the one obtained with the standard angle- 
average calculation. The amplitudes obtained here by retaining only the L = term 
were also observed to reproduce the spherically averaged results to excellent precision. 
This was used as a check of the computational algorithms. 

Values of the density-dependent G-matrix calculated with the spherical Pauli oper- 
ator were used for partial waves where 7 < J < 15, matrix elements that we needed to 
specify the long-range pion tail of the NN interaction. No partial waves with J > 15 
were included. The amplitudes of Eqs. (29)-(32) were then calculated at a number of 
values of and reproduced using a sum of Yukawa functions [5,14]. Because the L = 2 
contributions were much smaller than for L = 0, only L = 2 was considered and terms 
with L > 4 were ignored. 

The matrix inversion scheme described in Appendix A of Ref. |5] presumes that the 
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amplitudes associated with the Si 2 (q) and Si 2 (Q) are related according to G^f{ji~6) = 
— (— ) 5+T Gi^xiO) (as reported for E' and F' in Love and Franey For our present 

situation in which P is chosen to be parallel to (k+k'), this relationship is no longer valid 
(see the discussion in Appendix B). So we adopted an older fitting scheme in which the 
Gf ST (9) are calculated as a function of scattering angle and reproduced using standard 
least squares minimization techniques to determine the Yukawa expansion coefficients. 
Separate coefficients were obtained for Gj^{9) and Gj^(6), leading to the same quality 
of reproduction usually obtained for the spherical Pauli blocking case. 

4 Results for (p,p') reactions 

In the framework referred to by the name of Brueckner-Hartree-Fock (BHF), density- 
dependent effects on the interaction are included through Eq. (2), with the energy 
denominator properly modified by the presence of the nuclear matter mean field. Con- 
siderations of nucleons moving through nuclear matter as Dirac particles result in what 
is referred to as the Dirac-Brueckner-Hartree-Fock (DBHF) approach |2Sj- The latter is 
known to provide a realistic description of the nuclear matter equation of state |17j . 

It is useful to gauge the difference between the spherical and non-spherical Pauli 
operators against the background of these other contributions to medium modifications. 
To illustrate their effects, in Fig. 1 we show the cross section and analyzing power data 
for the 3~ state at 3.736 MeV in 40 Ca [21]. The beam energy for these measurements 
was 200 MeV. The short dashed curves contain no density dependence. Including only 
BHF effects produces the long-dashed curves; DBHF calculations produce the solid 
curves. As discussed in Ref. jS], the distortions are based on a folded optical potential 
that was constructed using the same effective interaction as the DWIA transition. The 
local nuclear matter density was unfolded from the longitudinal electron scattering form 
factor using the proton charge distribution. The structure formfactors replicate the 
inelastic electron scattering measurements for this same transition [21|. The calculations 
were made with the zero-range program LEA With the treatment of exchange 

discussed in Appendix B of Ref. [3] , the zero-range approximation provides an adequate 
description of natural parity transitions, as shown by the comparisons illustrated in 
Fig. 3 of Ref. [Sj. Using this method, we can maintain a high quality treatment of the 
transition formfactor. 

The most striking effects of the increasing repulsion in the nuclear medium lie beyond 
20° where they drive the analyzing power to more negative values and increase the cross 
section. With DBHF, very good agreement to the analyzing power data is obtained 
between 8° and 35° where the cross section is the largest and the DWIA should be at 
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e c.m. ( de 9) 

Stephenson, FIG. 1 

Figure 1: Measurements at 200 MeV of the cross section and analyzing power for the 
40 Ca(p,p') transition to the 3~ state at 3.736 MeV. The curves are DWIA calculations 
that include no medium effects (short dash), spherical BHF effects (long dash), and 
DBHF effects (solid). 
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its best. However, the cross section is overestimated by a factor of about 2, a problem 
of unknown origin with the DWIA that was circumvented for an empirical effective 
interaction by simply decreasing the normalization of the NN amplitudes at all densities 
[21]. If the calculated cross sections are lowered in this manner to agree with the data 
at the cross section peak, the DBHF curve provides the best reproduction of the shape 
of the angular distribution. 

In Fig. 2 we show for the same transition BHF-based predictions with non-spherical 
(solid) and spherical (dashed) Pauli operators. (At this time we do not have DBHF 
calculations with non-spherical Pauli blocking.) The differences are modest at all angles, 
being no more than a small fraction of the changes shown in Fig. 1. So given the quality 
of agreement with the DBHF calculations from Fig. 1, it is not possible to conclude 
whether the addition of non-spherical components to the treatment of Pauli blocking is 
required by the data. 

In Fig. 3, the non-spherical (solid) and spherical (dashed) BHF calculations are 
compared with measurements of the cross section [21], analyzing power [21], and spin 
rotation parameter Q(6) [2Z] for proton elastic scattering from 40 Ca at 200 MeV. To 
reduce complication, calculations without density dependence and DBHF calculations 
are not shown. Our earlier results |S] demonstrate that best agreement with elastic 
scattering is often obtained with just the BHF density dependence, a conclusion that 
has historically been verified many times [4]. Finite nucleus effects substantially reduce 
the size of the relativistic medium effects for elastic scattering [2H] • This reduction does 
not take place when the density dependence is obtained through calculations in infinite 
nuclear matter. Thus agreement between these BHF calculations and elastic scattering is 
good for both A(6) and Q(6) in the middle of the angular range of the data. Relativistic 
effects at angles below 15°, when included, do increase the analyzing power and shift the 
spin rotation parameter toward more positive values. Thus these differences where the 
elastic cross section is large can be reduced in the fuller treatment of the medium. The 
differences between the Pauli operators with and without non-spherical components are 
again too small to have an impact on the quality of the agreement with experimental 
data. 

Last, it is known that Pauli blocking effects become weaker as the bombarding energy 
goes up, finally disappearing near the pion production threshold. Figure 4 shows the 
non-spherical (solid) and spherical (dashed) BHF calculations for the 3~ state in 40 Ca at 
100 MeV. As expected, the differences are larger than in Fig. 2, but still remain smaller 
than the changes associated with the baseline medium effects or a DBHF treatment. 
(Earlier calculations j2H] that showed much larger effects at 100 MeV were in error.) 
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e c.m. ( de 9) 

Stephenson, FIG. 2 

Figure 2: Comparison of DWIA calculations with spherical (dashed) and non-spherical 
(solid) Pauli density dependence. The calculations at 200 MeV are for the 40 Ca(p,p') 
transition to the 3~ state at 3.736 MeV. 
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Figure 3: Measurements of the cross section, analyzing power, and spin rotation param- 
eter Q at 200 MeV for proton elastic scattering from 40 Ca. The folded optical model 
calculations are based on density-dependent interactions that include either spherical 
(dashed) or non-spherical (solid) Pauli effects. 
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Stephenson, FIG. 4 



Figure 4: Comparison of DWIA calculations with spherical (dashed) and non-spherical 
(solid) Pauli density dependence. The calculations at 100 MeV are for the 40 Ca(p,p') 
transition to the 3~ state at 3.736 MeV. 
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5 Conclusions 



A number of authors have examined the non-spherical Pauli blocking operator and cal- 
culated the changes to the G-matrix elements compared to the spherical approximation. 
We report here for the first time a formalism that allows the effects of the non-spherical 
operator to be included in a calculation of nucleon-induced elastic and inelastic scat- 
tering. We truncate the resulting G-matrix to remove the elements with J ^ J' or 
\i — £'\ > 2. This truncation involves terms that are small relative to the changes in 
the G-matrix elements with M, the projection quantum number for J. An expansion 
is introduced in terms of L = J — J' where the L = part essentially recovers the 
spherically- averaged Pauli blocking result. In this expansion we restrict our considera- 
tion to terms with L < 4. We wished to look in particular at polarization observables 
that might be sensitive to the M-dependence of the interaction matrix elements. 

Brueckner-Hartree-Fock calculations of the effective interaction including the non- 
spherical Pauli operator were made at 100 and 200 MeV. They were compared against 
other BHF-predictions obtained with the more common angle-average approximation, 
as well as DBHF calculations and those that contained no density dependence at all. 
The change from spherical to non-spherical in the treatment of the Pauli blocking op- 
erator produced changes that were small in comparison to the effects associated with 
typical BHF or DBHF medium modifications. Given the experimental errors and the 
likely larger uncertainties in the DWIA theory, the small size of the changes associated 
with adopting the non-spherical Pauli treatment cannot be deemed significant. These 
conclusions remained valid with changing bombarding energy. No polarization observ- 
ables (many were examined) were found to be particularly sensitive to the non-spherical 
treatment of the Pauli operator. Thus at the level at which present theory can provide 
an accurate model of nucleon-induced reactions, we find the spherically-averaged treat- 
ment of the Pauli projection operator adequate for intermediate-energy proton-nucleus 
scattering. 

The authors acknowledge financial support from the U.S. National Science Founda- 
tion under grant number PHY-0100348 (E.J.S.), the U.S. Department of Energy under 
grant number DE-FG02-03ER41270 (F.S.) and the U.K. Engineering and Physical Sci- 
ences Research Council under grant number GR/M82141. 
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A Symmetries of the G-matrix in infinite uniform 
nuclear matter 



We define the G-matrix as the solution of Eq. (2) and think of it as the set of matrix 
elements of an operator in the NN relative coordinate space for a fixed P and hp. 
This operator does not conserve the total angular momentum in the NN center-of- 
mass if P 7^ 0. However, the Pauli blocking operator Q is invariant under spatial 
reflections and time reversal and we will assume that the NN interaction V is also. 
These symmetries place important restrictions on the free-space G-matrix. Here we 
examine their consequences for the structure of the on-shell G-matrix in nuclear matter. 

If the two-body interaction V is invariant under time reversal, it is straight forward 
to show formally from Eq. (2) that the G(P) operator for fixed P satisfies 

KG(V)K~ l = G(P) f , (33) 

where K is the anti-unitary time reversal operator for the system. What this relationship 
means is that for any states \(f>) and \ip) 

(0 I G(P) I V) = m I G(P) I Ktf,), (34) 

where \Ktp) = K \ tp). 

Our angular momentum states are chosen to have phases such that 

K I (IS)JM) = (-l) J+M+e I (IS) J - M) . (35) 

The occurence of the phase (— l) e in Eq. (35) implies that we have not included i l factors 
with the orbital angular momentum wavefunctions Yg m - This is consistent with what is 
assumed in Eq. (10). 

Using the results of Eqs. (34) and (35) we obtain 



{(l'S)J'M' I G T (P) I (IS)JM) 

= ^y'+J'+M'+J+i+M^j _ M J G T (p) J ygjj, _ M>) 

In the case of no Pauli blocking, G conserves J and M and the G-matrix is indepen- 
dent of M. Together with parity conservation, this leads to the usual result that the 
G-matrix is symmetric in the angular momentum basis. 

When Eq. (36) is used in the definition of G LTS from Eq. (9), we find 

= (-D— ia^urAP), (37, 
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and using this in Eq. (25) (assuming k = k') gives 

Ggf s (k,k',P) = (-l) fc *Gg£(-k',-k,P) . (38) 

Eq. (38) is valid for arbitrary fixed P. 

For fixed P a consequence of parity conservation in Eq. (25) is 

Gjg£(k, k', P) = Gg£(-k, -k', P) , (39) 

Taken together, Eqs. (38) and (39) imply 

Gjg£(k, k', P) = (-l)^G^(k', k, P) . (40) 

Eqs. (38) and (39) are the generalization of the usual prescriptions for constructing G. 
They lead to the form given in Eq. (15) in the spherical Pauli-blocking approximation. 
For a general P, Eqs. (38) and (39) do not restrict G to have the form of Eq. (15) and 
there will be many other terms involving P, for example S^P). 

We can look at the consequences of Eqs. (38) and (39) in a different way. Eq. (40) is 
a relationship between two different sets of G-matrix elements corresponding to different 
ingoing and outgoing momenta with the same (arbitrary) P. We can change this into 
a relation between matrix elements with the same ingoing and outgoing momenta but 
with different P's. We do this by noting that, if k = k! interchanging k and k' (as in 
Eq. (40)) is the same as rotating them by 7r about k + k'. Eq. (40) is therefore equivalent 
to 

Gg£ (k, k', P) = (-l)**Gg£(tfk, Hk, HP') , (41) 
where R denotes the rotation by ir about k + k' and 

p' = ir x p . (42) 

We next use the fact that the G|J are components of an irreducible tensor of rank 
ks,qs- Therefore the condition Eq. (41), can be written [3U| 

Ggf s (k,k',P)) = (-l) fcs E^J4(k,k',P')P^ s (k + k',7r) , (43) 

i's 

where V> k q f q (k + k', re) is the rotation matrix corresponding to a rotation of rr about 
k + k'. 

A similar argument based on the parity condition, Eq. (39), and the fact that — k = 
R(tt, n)k and -k' = R(tc, n)k' gives jUl, 

Gf™(k, k', P)) = (-l) k * £ Ggg(k, k', P")V k l qs (n, tt) , (44) 

i's 
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where P" = Br 1 ^, n)P. 

In free space, or in the spherical Pauli blocking approximation, the G-matrix does 
not depend on the direction of P and Eqs. (43) and (44) become linear relationships 
between tensor components G%Jq of the same rank. Expressed in terms of components 
in the "standard" coordinate system, Eq. (44) (parity) gives 

Gg2 = (-l) ks+ * s Gt™ qs . (45) 

This simply states that for ks = 1 the vector G^ 3 is proportional to n and that the 
ks = 2 tensor is determined by 3 complex amplitudes instead of the general 5: 

r<LST _ n fiLST _ r<LST 

Lr 1n — U, wii — Lri i 

!° ' _ 11 _ 1-1 _ (46) 

riLST r* L ST _ _riLST riLST _ riLST \ > 

"20 ' "21 ~~ "2-1 ' "22 ~~ "2-2 ■ 

Hence for S = 1 a consequence of Eq. (45) in free space is that the G-matrix is deter- 
mined by 5 complex ^-dependent amplitudes. 

Equation (43) (or Eq. (38), time reversal) gives no extra information for ks = 1, but 
for ks = 2 gives the extra condition 



sin 0( ^3/2 G^ T - G^ T ) + 2 cos 9 G^ T = , (47) 

so that in fact there are only four independent G-matrix amplitudes. (That Eq. (47) 
is the only extra condition can easily be seen by expressing Eq. (44) in a coordinate 
system with z along k + k'. It is found that Eq. (44) implies G^f" 7 " = in this system 
|3()j.) This is the content of the form given in Eq. (15). 

For a general P in nuclear matter none of these restrictions will be true because 
Eqs. (43) and (44) relate tensor components with different P. However, the G-operator 
satisfies G(— P) = G(P) and hence if 

P" = ±P (48) 

and 

P' = ±P , (49) 

where P' and P" are defined in Eq. (42) and following Eq. (44), respectively, then 
Eq. (44) does relate tensor components with the same P and Eqs. (46) and (47), will 
be satisfied. 

This situation occurs when P lies in the scattering plane; both Eqs. (48) and (49) 
are satisfied if P is parallel to k + k'. This is the choice we make in this paper. With this 
choice the conditions, Eqs. (46) and (47), are satisfied and, equivalently, the G-matrix 
has the form of Eq. (15) on shell. 
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B Direct and exchange tensor amplitudes for arbi- 
trary P 

We can ensure that the amplitudes Gf ST (9) of Eq. (15) are anti-symmetrized by sub- 
tracting in each case the same amplitude with the additional phase factor (— ) 5+T and 
with k replaced by — k. This requires that any value of I that appears in a partial wave 
decomposition of G\ ST {ff) satisfies (— y +s+T = —1. 

The direction of P, the total momentum of the colliding nucleons defined in Eq. (1), 
remains the same for both terms in the anti-symmetrization subtraction. If we consider 
just the amplitudes associated with the tensor spin operators Si 2 (q) and Si 2 (Q), then 
in terms of the unsymmetrized amplitudes Gf ST (8): 



(k'|G(P)|k) = Gff s D T (cos9, Pq, PQ, Pn) S 12 (q)+G^ S X T (cos 9, Pq, PQ, Pn) S 12 (Q) (50) 



where the dot products indicate any scalar generated with the vector P. The amplitude 
to be subtracted is 



Note that the definition of the scattering angle becomes the complement of the original 
value. With the reversal of the direction of k, the vectors q and Q have interchanged 



G^(6) = G^ s J{cos6, P q, P Q, Pn) - (-) S+T G^ s /(- cos8, P Q, P q, -P n) (52) 



(k'|G(P)| -k) = G^ r (-cos#,P-Q,P-q,-P-n) Sl2 (Q) 
+ G£! T (-cos0,P-Q,P-q,-P-n) S 12 (q) . 



(51) 



places and magnitudes while P remains fixed. Thus, the direction of n = q x Q reverses. 
By combining terms, the anti-symmetrized form Gf ST (6) is given by 



for the coefficient of Si 2 (q) and 



G L A T (e) 




(53) 




(54) 
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This permits the same Yukawa expansion coefficients to be used for G^ s ^{6) and G^x" '($) 
provided that Q is used in place of q in the Yukawa expansion of the TX term ^3] • For 
the non-spherical calculations used in this paper, P is along the direction of Q, thus the 
dot products involving Q in Eqs. (52) and (53) remain. Due to the presence of different 
arguments in the TD and TX amplitudes, Eq. (54) is not necessarily satisfied. This is 
the case for the calculations reported here. 
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